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The two highest precision determinations of a s (M|), that based on the analysis of 
short-distance-sensitive lattice observables, and that based on an analysis of hadronic 
t decay data, have, until very recently, given results which are not in good agreement. 
I review new versions of these analyses which bring the two determinations into 
excellent agreement, and discuss prospects for additional future improvements. 
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I. INTRODUCTION 

Until recently, the determination of a s (M|) from perturbative analyses of short- 
distance-sensitive lattice observables (yielding 0.1170(12) [l|), and that from finite energy 
sum rule (FESR) analyses of hadronic r decay data (yielding 0.1212(11) j2]), produced 
central values which, though nominally precise, differed from one another by ~ 3a. In 
\0 ' the past year, this discrepancy has been removed by new versions of both analyses. In 

O 



what follows, I outline the important features of these updates which are responsible for 



this change. 



II. THE LATTICE DETERMINATION 



The original lattice determination |l) involved the perturbative analysis of various 
lattice observables, O k , computed using the MILC rif = 2 + 1 configurations. The D = 
expansions, computed to 3-loops for the MILC action by Mason et al. [l|, 0], have the 
form 

Ok = ^N<4(Qk) = D k a v (Q k ) <$i<Xv(Qk) (1) 

N=l M=0 

where ay is defined in the recent update ^ (HPQCD08), and Qk = d k /a is the relevant 
BLM scale. The c^\ z (equivalently, D k ,c[ \ <4 ), and d k are tabulated in HPQCD08. 
Regarding possible higher D contributions, (i) m^-dependent contributions were removed 
by extrapolation, using data, and (ii) non-perturbative (NP) m g -independent higher D 
contributions were initially assumed to be dominated by D = 4 gluon condensate terms, 
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which were then fitted and removed independently for each Ok- For the lattice spacings, 
a ~ 0.18, 0.12, and 0.09 fm, of the original analysis, the observed scale-dependence of 
the Ok could be reproduced only by also fitting at least one additional coefficient in 
Eq. ©13. 

The updated HPQCD08 [3] and CSSM [4j analyses incorporate data from new MILC 
ensembles with a ~ 0.15 and 0.06 fm, with one very new a ~ 0.045 fm ensemble also 
employed in HPQCD08. Useful cross-checks are also provided by differences in (i) the 
choice of coupling employed and (ii) the treatment of m q - independent NP contribution 
in the two re-analyses. The choice of coupling in HPQCD08 leaves residual perturbative 
uncertainties in the conversion from the V to MS scheme, that in CSSM in the effects 
of the truncated (3 function, which can be suppressed by focussing on finer lattices Q. 
HPQCD08 performs an improved treatment of m g -independent NP contributions, fitting 
a range of D > 4 forms to data, while CSSM restricts its attention to observables 
where the corresponding D = 4 contributions, estimated using charmonium sum-rule 
input for (a s G 2 ) jBf, are found to be small. Even with the finer lattice scales of the 
new analyses, it turns out that at least one additional coefficient in Eq. (TjQ) must be fit. 
The a s which results provides an excellent representation of the scale dependence of all 
the Ok employed. The results of the two re- analyses are in good agreement, and differ 
by only ~ lcr from the results of the earlier lattice analysis. The equivalent n/ = 5 
results, at the scale fi 2 = M§, are shown in Table H] for the three most- and four least- 
perturbative of the Ok studied in HPQCD08. Wu is the k x I Wilson loop and uq = W\l^ . 
Also shown is the corresponding quantity <5d=4, equal to the percent shift in the scale 
dependence between a ~ 0.12 and a ~ 0.06 fm resulting from first computing Ok using 
raw simulation values and then recomputing it after subtracting the known leading order 
m g -independent D = 4 contributions, estimated using charmonium sum rule input for 
{a s G 2 ) [5|. Sd=4 provides a measure of the expected importance of the m^-independent 
NP subtractions relative to the D = contribution from which a s is to be determined. 
We see that residual NP effects in Ok = log(Wn) are expected to be tiny; the resulting 
subtraction in fact produces a shift of only 0.0001 in a s {M\ ) 0]. The fact that the results 
for a s obtained after fitting and subtracting what are expected to be rather sizeable NP 
contributions to the most non-perturbative of the Ok agree so well with those obtained 
from analyses of the data for those Ok where these corrections are expected to be very 
small gives increased confidence in the treatment of such NP contributions and, even 
more important, enhanced confidence in the reliability of the results obtained from the 
most UV-sensitive of the Ok considered. 



III. THE HADRONIC r DECAY DETERMINATION 

In the SM, with Py/^. ud the / = 1 V or A current-induced width for r to hadrons, 
T e the r electronic width, y T = s/m 2 , Sew a known short- distance EW correction, and 
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TABLE I: a s (Mt) and 5d=a for the 3 least- and 4 most-non-perturbative of the 



O k 



a s (M|) a s (Ml) S D=4 
HPQCD08 CSSM 



logWi) 
log(W 12 ) 

log 



0.1185(8) 0.1190(11) 0.7% 
0.1185(8) 0.1191(11) 2.0% 

0.1183(7) 0.1191(11) 5.2% 



U,, ( ^) 0.1185(9) 
Zog ( -^f ) 0.1176(9) 



0.1171(11) 



N/A 
N/A 
N/A 
N/A 



32% 
53% 
79% 
92% 



R 



-phad 

V/A;ud — J- V/A;ud 



/r e , one has [8] 
dRv/A-,ud/dy T = I2n 2 Sew \V u d\ 2 w 00 (y T )py^ ] (s) 

-WLiy^Py/A^ 



(2) 



whore Pyj A {s) are the spectral functions of the spin J scalar correlators, Hy) A (s), of the 



1, V/A current-current 2-point functions, w 00 (y) 



-y) 2 (l+2y),w L (y)=2y(l-y) 2 
■ 2f 2 S(s — ml), making 

P v /' A '{s) accessible from the experimental distributions dRy/A ]U d/dy T (1,0]. For any s 
and any analytic w(s), the related finite energy sum rule (FESR) 



and, up to O [(m d ± m u ) 2 ] corrections, Py\s) 



and Pa\ s ) 



*0 



w ( s )Pv/~A l 



s)ds 



2rri 



w(s)U^\s)ds, 



(3) 



is satisfied. For large enough s 0) the OPE can be employed on the RHS. For typical 
w(s), and s above ~ 2 GeV 2 , the OPE is strongly D = dominated and hence largely 
determined by a s . The 5-loop version of the D = OPE series j9[ is employed in all 
recent r-based a s analyses 0, [9 14 1 . Use of polynomial weights, w(y), with y = s/s , 
helps to quantify higher D contributions, most of which must be fit to data, since (i) up 
to corrections of O {a 2 )-, the integrated OPE series terminates at D = 2N + 2 (with A" 
the degree of w(y)), and (ii) integrated OPE contributions of different D scale differently 
with s (D = 2k + 2 terms scaling as 

r decay determinations of a s have conventionally been based on a combined analysis 
of the s — m 2 km = 00, 10, 11, 12, 13, Wk m (y) = u>oo(j/) (1 — y) k y m "spectral weight 
FESRs" @, @, 0]. This analysis relies crucially on the assumption that D = 10, • • • , 16 
contributions, each in principle present for one or more of the Wkm employed, can, in 
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all cases,be safely neglected. A number of recent analyses employ either this strategy 
directly [2| or make use of nominal D = 6, 8 NP contributions extracted by doing so 0-11, 
141 ] . This assumption is, however, potentially dangerous since (i) a ~ 1% determination 
of a s (M|) requires control of D > 4 NP contributions to < 0.5% of the leading D = 
term and (ii) the fcm = 11, 12 and 13 FESRs have strongly suppressed D = OPE 
contributions. The validity of the assumption that all D > 8 contributions can be 
neglected was tested in MY08 [l2[ by (i) studying, as a function of s , the match between 
the variously weighted OPE integrals, evaluated using fitted OPE parameters, and the 
corresponding experimental spectral integrals, and (ii) using the same data and fitted 
OPE parameters as input to FESRs for different w(y) involving the same set of OPE 
parameters. The quality of the match produced by the results of the optimized Wk m 
analysis of the ALEPH data was found to be typically rather poor in the window ~ 
2 GeV 2 < s < ml, not just for the Wkm employed in the analysis, but also for other 
degree < 3 weights, whose OPE integrals should depend only on the D = 0, 4, 6, 8 OPE 
parameters included in the ALEPH fit [l2|. Similar discrepancies exist for the Wkm 
analysis based on the OPAL data, and for the BJ08 [llj treatment of the W np F ESR, 



which employs a different set of assumed values for the D = 6, 8 contributions [12 . 

In view of these problems, MY08 performed analyses based on alternate weights, 
w n{v) = 1 — W~[y + n~iU N i designed to suppress D > 4 relative to the leading D = 
contributions and hence optimize the determination of a s . (In terms of its size relative to 
the leading, a s -dependent integrated D = series, neglect of D > 8 contributions would, 
in fact, be between 7 and 814 times safer for the wn analyses than it is for the higher w^m 
FESRs of the conventional analyses [l2|. The fact that only a single D > 4 contribution, 
with D = 2N + 2, enters the wn FESR also simplifies fitting the corresponding OPE 
contribution.) One finds (i) excellent consistency among the a s values obtained using 
different W]y(y) and/or different channels (V, A or V+A); (ii) (as intended) a significantly 
reduced impact of D > 4 OPE contributions; and (iii) in contrast to the results of the 
combined {wk m } analysis, an excellent quality match between OPE and corresponding 
spectral integrals for other degree < 3 weights (including the kinematic weight woo) 
over the whole of the so window noted above. The results of MY08 are the only ones 
in the literature to satisfy this set of self-consistency constraints. Since (i) the {wk m } 
analyses, which should produce results in agreement with those of the corresponding 
wn analyses when using the same data, instead produce significantly larger a s , and 
(ii) (as shown in the Figures of MY08) the results of the {wk m } analyses, considered 
at lower sq, produce optimized OPE integrals not in agreement within errors with the 
corresponding experimental spectral integrals, and, moreover, significantly inferior to 
the matches obtained using the {wn} analysis fit parameters, we take the results for a s 
to be those obtained from the {wn} analysis of MY08. The central value, a s (m 2 T ) = 
0.321(5)(12) (where the first error is experimental and the second theoretical) is based 
on the CIPT treatment of the D = series, which yields better consistency among the 



results of the different wn FESRs than does the corresponding FOPT treatment [15 
This corresponds to 

a,(M|) = 0.1187(6)(15)(3), (4) 
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where the errors are, respectively, experimental, theoretical, and that due to evolution. 
The t result is now in excellent agreement with the lattice determination. Note that 
effects associated with the breakdown of the integrated OPE representation, if any, are 
not reflected in the errors quoted here. There are, in fact, no signs for the presence of 
such effects in the match between the optimized OPE and spectral integrals at present. 
Further discussion of this issue may be found in the next section. 

IV. DISCUSSION AND PROSPECTS 

Given the sizeable estimated D = 4 NP subtractions for the most non-perturbative 
of the lattice observables Ok, and the necessity of fitting analogous D > 4 contributions 
to data in an approximate way, the results corresponding to those Ok expected to have 
very small m 9 -independent NP contributions should be the most reliable sources of a s 
in the lattice analysis. The difference between the result obtained from the most UV- 
sensitive of the observables, Ok = log{Wu), a s (M^) = 0.1185(8), and that representing 
the average over the results for all the O k considered in HPQCD08, a fl (Mf ) = 0.1183(8), 
is however, in fact, very small on the scale of the uncertainties of the analysis. We note 
also that one would have to reduce the lattice spacing by roughly an order of magnitude 
before perturbative coefficients beyond 3-loops, which must be fitted at present, could 
be plausibly neglected. This is certainly not feasible. Fortunately, the fitting procedures 
employed, when applied to input pseudo-data generated using a known high-order input 
expansion, return the input a s value with very good accuracy. As a result, it appears 
highly unlikely that the errors associated with the fitting of the higher order coefficients 
play any significant role. The uncertainty in a s for what we would argue is the most fa- 
vorable of the Ok analyses, that based on Ok = log(Wu), as can be seen from the results 
quoted by HPQCD08, will thus be dominated by the uncertainties in the determination 
of the lattice spacings in physical units for the various ensembles, with comparable con- 
tributions coming from the uncertainties in r\ja and T\ itself jjjj]. It appears unlikely that 
these can be significantly reduced in the near future, so one should expect only minor 
improvements in the lattice determination, such as those that should result from access 
to a larger set of a ~ 0.045 fm ensembles. 

The current errors on the r determination are larger than for the lattice determination, 
but may be amenable to more significant near-term reduction. The major components 
of the quoted 0.012 uncertainty on a s (m*) are: 0.0084 from the maximum difference 
between CIPT and FOPT determinations over the W2, ■ ■ ■ ,wq set considered in MY08; 
0.0059 from the uncertainty in the 0.009(7) GeV 4 charmonium sum rule input employed 
for (a s G 2 /7r); and 0.0056 from the assumed 100% uncertainty on the FAC estimate for 
the 6-loop coefficient of the D = Adler function [||. We discuss these contributions 
briefly in what follows. 

BJ08 have raised an interesting question about the relative reliability of the FOPT 
and CIPT prescriptions for evaluating the truncated D = series. Based on a model 
reflecting known general properties of the divergent D = series, they argue that, de- 
spite better observed convergence behavior for the truncated CIPT series, the truncated 
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FOPT series might, nonetheless, better approximate the true result. Currently, it is 
known that, contrary to what one would expect from this model, a combined fit to a col- 
lection of degree < 3 weights using the truncated fifth order CIPT prescription yields a 
good combined OPE-spectral integral match, while the analogous combined fit using the 
truncated FOPT prescription does not [15|]. The analogous calculation, however, has yet 
to be carried out for the fully resummed BJ08 model. This work is in progress. An inter- 
esting additional observation is that the W2 case yields a much better agreement between 
the CIPT and FOPT treatments of the integrated D = OPE series than is found for 
other weights, and that this agreement persists over the whole of the 2 GeV 2 < Sq < m 2 . 
window. A much better than usual agreement between the FOPT and CIPT versions of 
the unweighted integrated truncated D = OPE series in the vicinity of the smallest 
terms is also seen in the BJ08 model, suggesting that the W2 FESR may be especially well 
suited to the determination of a s . Further work, however, is needed to have any realis- 
tic hope of significantly reducing the existing FOPT vs. CIPT prescript ion- dependence 
of the results for a s . The prescription-dependence uncertainty is, of course, intimately 
related to that associated with the uncertainty on the estimated 6-loop Adler function 
coefficient. 

Reduction in the error associated with the uncertainty on the input for (a s G 2 ) would 
require an improved determination of this condensate. Given the associated renormalon 
ambiguity, such a determination should be performed using the same hadronic r decay 
data, in an analysis with a consistently truncated D = series. Preliminary studies 
indicate that it is, in fact, likely that such an improved determination of (a s G 2 ) is 
possible. This issue turns out to be closely linked to the question of the level of residual 
duality violation (integrated OPE breakdown) present in the FESRs employed in the 
analysis, and the modelling of duality violating contributions to the physical spectral 
functions. 

A final potential error source for the r decay determination is residual integrated 
duality violation. While (i) it is known empirically that, at the scales employed in the 
r decay analyses, OPE breakdown is localized to the vicinity of the timelike point on 
the OPE contour [l6j], and (ii) the weights employed in the analysis described above 
all have a double zero at the timelike point, a physically well-motivated model for such 



duality violations [lOj], when fitted to the observed V and A spectral functions, allows 
duality-violating-induced shifts to a s (m 2 ) in the range 0.003 — 0.010. An additional 
uncertainty at the upper end of this range would not be negligible on the scale of the 
other errors quoted above. Preliminary investigations, however, indicate that significant 
further constraints can be placed on this model and that, when they are, (i) the allowed 
duality-violating-induced shifts lie at the low end of the range quoted above, and (ii) 
an improved simultaneous determination of the gluon condensate is almost certainly 
possible. Further work is required before more quantitive statements on these issues can 
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be made. 
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